Extended scaling for ferromagnetic Ising models with zero-temperature transitions 
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We study the second-moment correlation length and the reduced susceptibility of two ferromag- 
netic Ising models with zero-temperature ordering. By introducing a scaling variable motivated by 
high-temperature series expansions, we are able to scale data for the one-dimensional Ising ferromag- 
net rigorously over the entire temperature range. Analogous scaling expressions are then applied to 
the two-dimensional fully frustrated Villain model where excellent finite-size scaling over the entire 
temperature range is achieved. Thus we broaden the applicability of the extended scaling method 
to Ising systems having a zero-temperature critical point. 
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Studying the critical behavior of systems that order 
at zero temperature is challenging because the typi- 
cally used Monte Carlo methods are generally unable 
to probe the critical behavior close enough to the zero- 
temperature critical point for traditional 1 finite-size scal- 
ing approaches to yield precise critical parameters. It is 
thus necessary to either incorporate scaling corrections 2 
or find better approaches to scale the d&t&££ 

Using the intuition gained from high-temperature se- 
ries expansions (HTSEs) a scaling approach has been in- 
troduced with the aim of extending the validity of crit- 
ical scaling expressions to temperatures well above the 
critical region. So far, the approach has been applied 
to a number of model systems having finite ordering 
temperatures^^ Inherently, this approach is ideal to 
study systems which order at zero temperature — such as 
a magnetic system below the lower critical dimension — 
since there only temperatures above the critical point can 
be accessed numerically. Thus, an important first step in 
adapting the extended scaling approach to these systems 
involves the appropriate choice of a scaling variable. 

In this paper we derive extended scaling relations for 
two sample ferromagnetic models with no disorder in the 
interactions and which only order at zero temperature. 
First, we study the exactly solvable one-dimensional (ID) 
Ising fcrromagnct and then use insights obtained to ana- 
lyze the nontrivial two-dimensional (2D) fully frustrated 
Villain Ising models Data generated using Monte Carlo 
simulations for large system sizes and very low tempera- 
tures validate our scaling approach. 

The paper is structured as follows. In Sec. [IT] we dis- 
cuss the extended scaling approach and how to adapt 
it to systems ordering at zero temperature, illustrating 
the results with the one-dimensional case in Sec. IIIII In 
Sec. IIVI we introduce the (Villain) fully frustrated Ising 
model and present details of our numerical calculations. 
The numerical data are then analyzed with the extended 
scaling approach, followed by concluding remarks. 



Conventionally, at a continuous phase transition the 
power-law critical behavior of any observable O in the 
thermodynamic limit can be written asi 



0[t(T)} ~ f 



(1) 



with t = (T — T c )/T c as a scaling variable, T the tem- 
perature and T c the critical temperature at which the 
phase transition occurs. The exponent y describes the 
"strength" of the divergence at T c . Alternatively, other 
critical variables which yield the same limiting behavior 
at criticality such as 



r=(T-T c )/T=[l-/3//? e ], 



(2) 



where [3 = 1/T, can be used. For certain models further 
scaling variables have been introduced, e.g., 



r s = i[sinh- 1 (2/3)-sinh(2/3)] 



for the 2D Ising ferromagnet 



10.11.12 



= [1 - (P/Pcf 



(3) 



(4) 



for spin glasses with zero-mean symmetric interaction 
distributions All these scaling variables are propor- 
tional to (T - T c ) for T -» T c . However, T is not a 
sensible scaling variable for a zero-temperature (T c = 0) 
transition^ Cardy et al*£ gave a renormalization-group 
technique rule for low-temperature limiting scaling vari- 
ables at zero-temperature transitions. For example, for 
the Potts model studied by Cardy et al. the appropriate 
renormalization-group "temperature" scaling variable is 
proportional to exp(— (3). By analogy, for a ID ferromag- 
netic Ising model for which T c — the scaling variable 
should be proportional to cxp(— 2(3)*^ 

Scaling expressions may also include temperature- 
dependent prcfactors, which are noncritical but can be 
relevant in analyses that include a range of temperatures 
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far above (or below) T c . For instance, the reduced sus- 
ceptibility x(P) measured in numerical simulations (see 
below) is related to the thermodynamic susceptibility 
Xth(P) (which is the physically measurable observable) 
through x(P) = Xth{P)/P- Thus, critical behavior of the 
form x{P) ~ i~ 7 implies Xth(/3) ~ Pt -7 , i.e., with a pref- 
actor 0. 

The two thermodynamic limit obscrvables that we dis- 
cuss here are the ferromagnetic reduced susceptibility \ 
and the second-moment correlation length £. The ferro- 
magnetic reduced susceptibility is given by 



where 



X = N(m 2 ), 



1 N 



(5) 



(6) 



is the magnetization per spin, ./V is the number of spins in 
the system and (• • ■ ) represents a thermal average. The 
second-moment correlation length is given by 



M2 

zx 



1/2 



where 



N 

1-12 = r ij( S i S o) 
i,j=l 



(7) 



(8) 



is the second moment of the correlation function, z is the 
number of nearest neighbors and is the distance be- 
tween spins i and j J£ For a hypercubic lattice z = 2D, 
where D is the space dimension. Note that numerically 
we measure the finite-size correlation length [which is 
equivalent to the expression presented in Eq. |7|)] as 



2sin(|k m |/2) 



X 



X(km) 



- 1 



1/2 



(9) 



where k m = (27r/L,0) is the smallest nonzero wave vec- 
tor (here in 2D), and x(k) is the wave-vector-dependent 
reduced susceptibility: 



1 N 
(k) = _ £ {SiSj) 



(10) 



HTSEs of the Ising ferromagnct in large space dimen- 
sions (i.e., in the mean-field regime )i show that simple 
relations for the reduced susceptibility, namely 



x(P) 



(11) 



and for the second-moment correlation length defined in 
Eq. 



are exact for all T > T c = 1. Thus, in this limit with the 
scaling variable t [Eq. (J2J], the critical power laws for the 
reduced observables Xth (P) I P and £ (P) / P 1 ' 2 hold exactly 
over the entire range of /3 from /3 C to zero. In finite- 
dimensional ferromagnetic systems — if the same basic 
variables and expressions are used^&i with the modifica- 
tion necessary to give the right high-temperature limits — 
one obtains "extended scaling" equations in which the 
leading terms are 



and 



X {P) = C X T-T + (1 - C x ) 



m = P 1/a [Cer-' + {i-Ct)] 



(13) 



(14) 



In Eqs. (fT3")) and (fT4]) C x and Cj are critical amplitudes 
and 7 and v are the standard critical exponents^ With 
the appropriate critical parameters (critical temperature, 
critical exponents, as well as critical amplitudes) these 
expressions are exact by construction at the (3 — > f3 c and 
P — > limits. Elsewhere, the expressions are not ex- 
act but have been shown to give good approximations to 
the true behavior for the entire paramagnetic tempera- 
ture region. By introducing small correction terms these 
approximations can be improved considerably. 

The expression for finite p c given in Eqs. (fT3|) and (jT4j) 
cannot be used for systems with T c = because P c = oo. 
In Sec. |TTT] we present an "extended scaling" approach 
tailored to systems having T c = and a unique nearest- 
neighbor interaction strength | Jy | (in this case | Jy | = 1 
V i, j and no bond disorder)^ We first present simple 
exact expressions for the ID Ising ferromagnet for which 
the scaling variable 



T t {P) = l-tanh03) 



(15) 



works well. This is consistent with the Cardy et al. rule 3,5 
because r t is equal to 2exp(— 2/3) at low temperatures; 
but, like (1 — P/Pc), r t tends to 1 for T — > oo. In the light 
of this result we then apply the same approach using r t 
to the nontrivial 2D fully frustrated Villain model. Our 
analysis shows that the extended scaling scenario with 
r t as a scaling variable gives an excellent account of the 
behavior of the correlation length and reduced suscepti- 
bility extrapolated to infinite size over the entire temper- 
ature range. 



III. ONE-DIMENSIONAL ISING MODEL 

To motivate the scaling expressions for ferromagnetic 
Ising models with zero transition temperature, we use as 
a toy model the one-dimensional Ising ferromagnet 



H 



ID 



(16) 



(12) 



with Ji^+i = 1 for all nearest neighbors i and i + The 
model orders onljsii at T — and expressions for £(/?) 
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FIG. 1: (Color online) Scaled ferromagnetic reduced suscep- 
tibility [Eq. ((5}] for the ID Ising ferromagnet according to 
Eq. p5|l . The data scale perfectly and thus validate the de- 
rived scaling expressions. The inset shows the unsealed data 
for different system sizes, as well as the thermodynamic limit 
(thick gray line). 



and x(/3) in the infinite-size limit are easily calculated 
from HTSE. The reduced susceptibility can be expanded 
as 

X(J3) = 1 + 2[tanh(/3) + tanh 2 (/3) + tanh 3 (/3) + . . .] (17) 

and the second moment of the correlation is 

fj, 2 (P) = 2[tanh(/3) + 2 2 tanh 2 (/3) + 3 2 tanh 3 (/3) + ...]. 

(18) 

The second-moment correlation length is then given by 
Eq. ([7]) with z = 2 in ID. Using the mathematical iden- 
tities 



X> n = T_^ and E ; 



71=1 



(1-x) 3 ' 



(19) 



the exact expressions for susceptibility and correlation 
length are thus 



X ((3) = exp(2/3) 



1 - tanh(/3) 



- 1 



(20) 



and 



m = i[exp(4/3)-l]V2 

tanh 1/2 (/3) 
1 - tanh(/3) 



(21) 



Note that these expressions are valid for the entire tem- 
perature range. 

Equations (|2"D|) and (f2"Tj) are of the extended scaling 
formi [see Eqs. (H3J) and (|T3|)] with r t [Eq. fT5])] replacing 



t in the extended scaling expressions for the ferromagncts 
with finite ordering temperatures. Finally, temperature- 
dependent effective exponents can be defined as 



7 (/3) = -dlog[ X ]/dlog[l - tanh(/3)] 



(22) 



and 



v((3) = -dlog[£(/3)/tanh 1/2 (/3)]/dlog[l - tanh(/3)] . 

(23) 

In the limit T — > T c = the critical exponents are thus 

7c = y c = 1- 

For the Ising ferromagnet in one space dimension with 
linear extent L = N the Fisher finite-size scaling rule^ 
for an observable 

0(L,P)~L»l v T[Llt{P)\ (24) 
when applied to the reduced susceptibility leads to 



L/tanh(py/ 2 



L[l - tanh(/?)] 



tanh 1/2 (/3) 



L 



(25) 



In Fig. [T] we illustrate the previously-derived scaling re- 
lations with data for the reduced susceptibility for finite 
system sizes. The data are obtained by starting with the 
partition function Z = + for a one-dimensional 
system of L spins in a field ifjM with 



A- 



cosh(/?ff ) ± \J cosh 2 {(3H) - 26-^ sinh(2/3) 



(26) 

To obtain the thermodynamic susceptibility Xth(/3), we 
perform a second-order derivative of the free energy per 
spin, F = — In Z(L), with respect to H, subse- 
quently setting H = 0. The raw data for x = Xth/P 
(inset) are scaled according to Eq. (|25p. The scaling is 
perfect. 



IV. TWO-DIMENSIONAL VILLAIN MODEL 

The two-dimensional fully frustrated Ising model, or 
Villain model j£ consists of Ising spins on a square lattice 
with nearest-neighbor bonds \Jij\ = 1; in the x direc- 
tion all bonds are ferromagnetic, while in the y direction 
columns of bonds are alternately ferromagnetic and an- 
tifcrromagnctic. The Hamiltonian is thus given by 



H 



JijSiSj 

(i,3) 



(27) 



where Si = ±1 represent Ising spins on a square lattice 
with N = L 2 spins. The system is fully frustrated; i.e., 
the product of the bonds around each plaquette of the 
system is negative, 



-i . 



(28) 



4 



TABLE I: Parameters of the simulations. L denotes the sys- 
tem size, N Ba is the number of independent runs to improve 
the statistics and iV sw is the total number of Monte Carlo 
sweeps performed in a single run for each of the 2Nt replicas. 
T m in and T max are the lowest and highest temperatures sim- 
ulated, respectively, and Nt is the number of temperatures 
used in the parallel tempering method. 
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30 
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96 


100 


2097152 


30 


0.1 


3.0 



The model does not order at a finite temperature^ but 
exhibits a critical point at zero temperature with diverg- 
ing ferromagnetic reduced susceptibility and a ground- 
state degeneracy which grows exponentially with the sys- 
tem size. 

For the scaling analysis we compute the reduced sus- 
ceptibility [Eq. ([5])] and the finite-size second-moment 
correlation length [Eq. ©]. The simulations are done 
using exchange (parallel tempering) Monte Carlo 20 ' 21 i 22 
and the simulation parameters are presented in Table [U 
Equilibration is tested by a logarithmic binning of the 
data. Once the last three bins for all observables agree 
within error bars the system is considered to be in ther- 
mal equilibrium. We use periodic boundary conditions 
to reduce finite-size corrections. 

Forgac*^ showed analytically that the limiting low- 
temperature thermodynamic behavior of the correlation 
length of the 2D fully frustrated Villain model is strictly 
exponential, i.e., £(/?) ~ exp(2/3). Furthermore, the crit- 
ical exponent rj describing the decay of the correlation 
at T c is exactly 1/2 such that in the low-temperature 
limit, using X (P) ~ £,{PV ,V [Eqs. JB) and and 
the standard scaling relation 7 = (2 — rf)v, we obtain 
x{0) ~ £(/3) 2 ~ v = exp(3/3). Based on an analysis of the 
size dependence of the energy by Lukic et al^ it has been 
conjectured that the low-temperature limit for the cor- 
relation length is exactly £(/3) = (1/2) exp(2/3). No full 
HTSE study seems to have been carried out to date; how- 
ever, by inspection, the leading HTSE terms for the re- 
duced susceptibility are x(/3) = l + 2f3+0(f3 2 ) and for the 
second moment of the correlation length /X2 = j3+ O(0 2 ). 

Scaling dimcnsionless ratios of finite-size data for the 
correlation length [£(2L, /?)/£(L, (3)] and susceptibility 
[x(2L, f3)/x(L, (3)] vs the two-point finite-size correlation 
length divided by the system size [£(L,@)/L], which is 
also a dimensionless quantity^ yields unique curves de- 
pending only on the universality class if there are no 
finite-size corrections to scaling. For the system sizes 
studied, the Villain model shows weak corrections to scal- 
ing. This can be seen in Figs. [5] and [3] where the ratios 
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FIG. 2: (Color online) Scaling ratio £(2L, /?)/£(£, P) for the 
2D Villain model as a function of exp[— L/£(L, f3)] for different 
system sizes L. The dashed horizontal line corresponds to the 
exact infinite size value £(2L, f} c )/£(L, f} c ) = 2 at the critical 
point. The vertical line corresponds to the estimated infinite- 
size limit i{Po)/L = 0.488 (see text). 



are shown as functions of exp[— L/£(L, /3)]£ In principle, 
it should be possible to use the ansatz of Calabrese et 
al.-r^ i.e., 



£(2£,/?) 



= F[L/aL,f3)]+L-»g[L/aL,P)}, (29) 



where uj is the finite-size scaling correction exponent and 
T and Q are scaling functions. (Similar relations ap- 
ply for the susceptibility \-) As shown below when 
T — > 0, £00 (/3) diverges until £00 (/3) 3> L; the values 
of the observables saturate at £(L, 0) — ► £(L, C = 00) 
and x(-^)/3) —* x(L,/3 c )- The lowest temperature at 
which the simulations have been carried out is T = 
0.1. At this temperature we find £(L = 00,/?) ~ 
2.8 ■ 10 8 , thus for all L studied £(L = 00, 0) > L 
and we can take the measured values of observables 
at all L as good approximations to the T = value. 
Hence, for (3 — > (3 C the pre-factor sin(|k m |/2) is the 
only L-dependent factor in Eq. ©, which leads to 
£(2L,(3 C )/£(L,(3 C ) — ► 2. Furthermore, because 77 = 1/2, 
X (2L,p c )/x(L,p c ) -» 2 2 -" = 2.82843... exactly. Fig- 
ure [3|a) shows x(2L,@c)/x(L,@ c ) and Fig. fj|b) shows 
£(2L, P C )/£(L, P c ) against 1/L; it can be seen that the 
ratio behaves approximately as x(2L, ( 5 c )/x(i,/3 c ) ~ 
2.8284 + 0.U/L [£(2Z, C )/£(L, /3 C ) ~ 2 - 0.24/L] show- 
ing that the correction exponent for the leading correc- 
tion at large L can be plausibly taken as ui ~ 1 with 
further terms appearing at smaller L. In panel (c) of 
Fig. g] we show data for X (L, Pc)/L 2 -' n = x{L, f3 c )/L 3/2 
against 1/L. Fitting the data assuming ui = 1 gives the 
large-size limit x{L, (3 C )/L 3 / 2 w 0.585(1) - 0.05/L. In 
a similar way we find the approximate limiting value 
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FIG. 3: (Color online) Scaling ratio x( 2L , P)/x( L , ft) as a 
function of exp[— L/£(L, /?)] for the 2D Villain model for dif- 
ferent system sizes L. The dashed horizontal line corresponds 
to the exact infinite-size value x(2£, Pc)/x{L, Pc) = 2 3,/2 
at the critical point since in general x{Tc,L) ~ L 2 ~ v and 
77 = 1/2. The vertical line corresponds to the estimated 
infinite-size limit £(/3 c )/L = 0.488 (see text). 



of £(L,/3 C )/L ks 0.488(1) + 0.1/L [panel (d)] and of 
the Binder cumulant g{L,f3 c ) ~ 0.691(1) + 0.10/L (not 
shown) . 

The analysis of the data at other temperatures is also 
consistent with a leading correction with an exponent 
uj *~ 1 plus further correction terms for smaller L. At 
all temperatures studied the difference between the es- 
timated infinite-size values for the observables and the 
measured large- L values are always less than 0.5% of the 
measured large- Z/ values. 

Inspired by the results on the ID Ising ferromagnet 
outlined in Sec. IIHI with r t = 1 — tanh(/3) as a scaling 
variable we now test an analogous scaling of the data for 
the 2D fully frustrated Ising model. 

The critical exponents for the 2D fully frustrated Ising 
model [Eqs. (22]) and (221)] are j c = 3/2, u c = 1, and 
7] c = 1/2 J£ We thus construct trial expressions for the 
different observables as follows: 



Xff(P) =C x [1- tanh(/3)]" 3 / 2 + (1 - C x ) 



and 
£ff(/3) 



tanh 1/2 (/3) 



1 - tanh(/3) 



(1-Cf 



(30) 



(31) 



where the critical amplitudes C x and are the only ad- 
justable parameters. It turns out that for the correlation 
length, the expression with = 1.00(1), i.e., 



£ FF (/3)/tanh 1/2 (/3) = [1 - tanh(/3)]~ 



(32) 



gives an excellent overall fit to ^ 00 (/3)/tanh 1 ^ 2 (/3), which 
is the normalized infinite-size limiting curve estimated 
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FIG. 4: (Color online) Scaling ratios x(2£, P)/x{L, P) [panel 
(a)] and£(2L,/3)/£(L,/3) [panel (b)] plotted against 1/L. The 
dashed line corresponds to x(?L, P)/x(.L,0) ~ 2 3/2 + 0.U/L 
(£(2L,P)/£(L,f3) ~ 2-0.24/L). The full symbol corresponds 
to the exact thermodynamic value 2 3 / 2 in panel (a) and 2 in 
panel (b). Panel (c): x(L,/3)/L 3/2 vs 1/L. The dashed line 
corresponds to x(L, (3)/ L 3/2 ~ 0.585 - 0.05/L. Deviations 
appear for smaller values of L. Panel (d): £(L,/3)/L plotted 
against 1/L. The dashed line corresponds to ^{L,[3)/L ~ 
0.488 + 0.1/L. The data thus suggest that a corrections to 
scaling exponent ui ~ 1 is plausible. All data are for T — 0.10. 
Note that the data point for L — 96 is generally a bit high 
possibly due to the small statistics used in the simulation. 
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FIG. 5: (Color online) Normalized correlation length 
C(/3)/tanh 1/2 ( / 3) for the 2D fully frustrated Ising model for 
different system sizes L as a function of 1 — tanh(/3). The thick 
line corresponds to the extended scaling correlation length ex- 
pression Eq. (|32p . 



from the scaling curves; see Fig. [SJ Over the entire tem- 
perature range the maximum difference between the fit 
and the numerical curve is approximately 0.5%. The 
expression in Eq. (|3"2")l with = 1 is identical to the 
exact expression for the ID fsing ferromagnet. In the 
low-temperature limit with C% = 1, £ff — > (1/2) exp(2/3) 
meaning that the present data and analysis are consistent 
with the Lukic et al. conjecture 2 ^ within numerical pre- 
cision. In the high-temperature limit £ff — ► ft 1 ^ 2 , which 
is consistent with the first term of the HTSE for £(/3). 

For the reduced susceptibility the fits to the numerical 
data for x°o(/3) with Eq. ([50]) indicate that C x is equal 
to 2.00(5); see Fig. [H] The fit in the higher-temperature 
range can be improved further by a correction term cho- 
sen so that there is an exact agreement between the high- 
temperature limit obtained from the first two terms in 
the HTSE, namely, X {P) = [1 + 2/3 2 + . . .] as /3 ^ 0. We 
thus obtain: 

Xff(/3) = 2.0[l-tanh(/3)]- 3/2 -2.0+[l-tanh(/3)]. (33) 

With the extended scaling expressions given above 
[Eqs. (|32|) and (|33|) ] the standard Fisher finite-size scaling 
(FSS) [Eq. ([M])] is modified (see Refs.iandifor details). 
For the finite-size correlation length we thus obtain from 
Eqs. §M§ and 



L[l - tanh(/3)] 
tanh 1/2 (/3) 



= V 



FF 



(34) 
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FIG. 6: (Color online) Data for the susceptibility x(^> P) °f 
the 2D fully frustrated Ising model for different system sizes 
L as a function of 1 — tanh(/3). The thick line corresponds to 
the extended scaling expression in Eq. (|33|l . 



tain from Eqs. (|3T?j) and (|52|> 

[X(L,P) + 2 - (1 - tanh(/?))] 





[L/tanhiP) 1 / 2 


3/2 




L[l - tanh(/3)] 


— ^ X 


tanh 1/2 (/3) 



FF 



L 



(35) 



A finite-size scaling analysis of the data for the second- 
moment correlation length and the susceptibility using 
Eqs. ([34"]) and ([35)1 is shown in Figs. [7] and [3 respectively. 

The scaling curves have a characteristic form. Quite 
generally, at small £,00/ L, £(L,(3)/L = £(L = oo,0)/L 
so the log-log plot of, e.g., Fig. [7] is initially a straight 
line of slope 1 passing through the point [1,1]. At the 
large £,oc/L limit (which is equivalent to T = for all 
L) the curves tend to plateau values, K x = x(L)/L 3 ^ 2 = 
0.585(1) and K% = £{L)/L = 0.488(1) estimated above. 
If we ignore the marginal corrections to finite-size scaling, 
the crossover can be expressed phenomcnologically as 



£(L,0)/L 



Koo (£)/£] 



1 i/h 



1 + (1/K 6 ri[^(p)/LY 



(36) 



where is a crossover exponent. In the present case 
?s 2.5. For the reduced susceptibility the initial small- 
£00/ L behavior is X (L,(3)/L 2 ^ ~ [£«,(/?) / X]( 2 "") and 
the analogous phenomenological crossover equation is 



X (i,/3)/L 2 -" = A 



1 + (A/K X M^(P)/L}^'V) 



(37) 



whereas for the normalized reduced susceptibility we ob- 



where z x is the crossover exponent, A is a constant, and 
K x is the plateau value. The phenomenological fit values 
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FIG. 7: (Color online) Finite-size scaling of the two-point 
correlation length data of the 2D fully frustrated Ising model 
using the extended scaling expression, Eqs. (|32[) and (|34ll . 
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FIG. 8: (Color online) Finite-size scaling of the susceptibil- 
ity of the 2D fully frustrated Ising model using the extended 
scaling expression for the normalized susceptibility Xn(L,f3), 
Eqs. (|33]l and ((35]) . 

in the present case are z x k, 2.0 and A « 1.95. For both 
observables the fits with crossover are of excellent quality. 

V. SUMMARY AND CONCLUSION 

We have presented scaling expressions motivated via 
high-temperature series expansions which extend the 
scaling functions across the whole temperature range for 
systems which order at zero temperature. 

For the ID Ising fcrromagnet we can derive exact ex- 



tended scaling expressions of the form = 2r t _1 - 1 
and £(/3) = t^ 1 tanh 1//2 (/?) with the scaling variable 
r t = 1 — tanh(/3), critical exponents "f c = v c = 1 [defined 
via Eqs. (|22|) and ([23]) ]. and critical amplitudes C x = 2 
and Cf = 1. 

From the insights gained from the study of the ID fer- 
romagnet, we use the same temperature variable Tt to 
analyze numerical data for the 2D fully frustrated (Vil- 
lain) model. The exact critical exponents are knownJ^ 
7 C = 3/2 and v c = 1. We find that for the second- 
moment correlation length £(/3) = r t _1 tanh 1 / 2 (/?) with 
C*£ = 1 just as for the ID ferromagnet. Furthermore, this 
result is consistent within numerical accuracy with the 
low-temperature-limit conjecture of Lukic et al^ that 
£(/?) -> (l/2)cxp(2/3) for T -> 0; however the present 
expression covers the entire temperature range. The ap- 
proximate expression for the susceptibility of the Villain 
model [Eq. ([33])] with critical amplitude C x = 2.00(5) 
is in good agreement with the numerical data over the 
entire temperature range covered. 

Summarizing, the temperature dependence of observ- 
ables above a ferromagnetic transition (including T c = 
transitions) can be written in terms of generic extended 
scaling forms£&£ expressed to leading order as 

ttx)=x- 1 / 2 [Cs(l-x)- v + (l-C s )} (38) 

and 

x(x)=C x (l-x)-^ + {l-C x ) (39) 

with the scaling variable x and critical parameters de- 
pending on the system studied. The expressions are ex- 
act by construction in the limits (3 — » (3 C and (3 — * if the 
critical parameters are known. For a ferromagnet with 
T c > x = j3 1 (3 C . Note that in the high-dimensional limit 
Eqs. (]38]l and (]39]l are exact.— In finite dimensions (but 
with nonzero T c ) the expressions remain as good approx- 
imations over the entire temperature range. For the two 
ferromagnets with = 1, T c = 0, and no bond dis- 
order, x = tanh(/3) and thus rt = [1 — tanh(/3)] replaces 
[1 — (3/(3 c ] as the scaling variable. Effective exponents are 
defined through Eqs. (f^]) and P3"]h These relations are 
validated with numerical data on the 2D fully frustrated 
Ising model. 

There are numerous possible candidate systems to 
which this approach should in principle be applicable 
mutatis mutandis. These include for instance the fam- 
ily of fully frustrated 2D systems studied by Forgacs^ 
the 2D three-state Potts antifcrromagnet^i^, the 2D 
Ising antiferromagnet on a triangular lattice^ the 2D 
a modelsj^i as well as 2D Heiscnbcrg models<21 An 
interesting further step would be to determine scaling 
expressions for the two-dimensional bimodal Ising spin 
glass with | Jij | = 1 but with random signs for the in- 
teractions, which also orders at zero temperature. In 
that case the critical behavior of the model is highly 
controversia l 28 ' 29 ' 30 i 31 ' 32 i 33 and current data at finite 
temperature do not have the necessary quality within 
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"traditional" scaling approaches to determine the true 
nature of the transition. 
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